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Abstract 



o 

Q 

We present a study of the decay of metastable states of a scalar field via ther- 

■ mal activation, in the presence of a finite density of fermions. The process we 

■ consider is the nucleation of "drop/ets" of true vacuum inside the false one. 

(N 



We analyze a one-dimensional system of interacting bosons and fermions, con- 
sidering the latter at finite temperature and with a given chemical potential. 
As a consequence of a non-equilibrium formalism previously developed, we 



D 

,1^ • obtain time-dependent decay rates. 

> 

X 

^ . I. INTRODUCTION 

The dynamics of metastable systems is a subject that may be found in all realms of 
Physics. In Condensed Matter Physics, in polymer transitions and liquid mixtures P,^; 
in High Energy Physics, in the problem of baryon number violation [^] and in the formation 
of chiral condensates in heavy ions collisions |^ ; in Astrophysics, in "nucleating planets" ; in 
Cosmology, in the well known Inflationary Models |]^. All these examples have at least one 
thing in common: the occurrence of metastability, which requires non-equilibrium methods. 

As a model for learning something about the qualitative features of the decays of 
metastable states, we will study a one- dimensional system of interacting fermions and bosons 



1 



that starts in a metastable vacuum and gradually decays to the true one. The process con- 
sidered here will be the nucleation of "bubbles" of true vacuum inside the false one via 
thermal activation. Our main purpose will be the analysis of the stability of these "bub- 
bles" and the calculation of the decay rate as a function of time, in the presence of a finite 
density of fermions at finite temperature. As already shown in reference the decay rate 
is time-dependent as opposed to the time-independent results previously found in the liter- 
ature This is a direct consequence of the use of the non-equilibrium formalism. The 
results presented here appear to be in qualitative agreement with some experimental results. 

The introduction of fermions that interact with the scalar field has two remarkable fea- 
tures. The first, already used in reference is the preservation of the functional form of 
the bosonic static solution of the equation of motion for the case with no fermions. Through 
the inverse scattering methods used here, this fact is rigorously demonstrated. The second 
is the appearance of metastability within metastability. Besides the metastability which 
is inherent to our choice of the metastable minimum of the effective bosonic potential, we 
find that the extrema of such a potential may be metastable "bubbles" . Thus, throughout 
this article, the metastable state (for example, the first excited state of a cis-polymer) may 
develop defects which correspond to either unstable or metastable "bubbles" . 

The paper is organized as follows: in section |T|, we briefly review the results for just 



a scalar field, which already appeared in the literature; in section 111, we include fermion 



fields and analyze their influence on the effective theory for the bosons; in section we 
add the effects of flnite temperature and of a chemical potential for the fermions; in section 
0, we present some comments about the method and our conclusions. Some mathematical 
definitions are left for a final Appendix. 

II. THE BOSONIC CASE: A SUMMARY 

In this section we summarize results concerning the decay of metastable bosonic systems. 
Most of the material presented here may be found in various recent papers MM- 
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A. The Lagrangian 



Our bosonic system is represented by a real one-dimensional scalar field, The 
dynamics is given by a Lagrangian with the following form 

>C = - A(0 - 0_) V(0 - 0*) (2.1) 

where 

The Lagrangian (|2.1| ) has an asymmetric double well potential,\^(0) = X{(f) — </)_)^0(0 — 
(f)*), with V"{(f)) = m? , where im? is the squared mass of small oscillations around the local 
minimum 0_ (see Figure 1) and A is a coupling constant. 



B. The "sphaleron" solution 



The equation of motion associated with ( |2.1| ) has a static solution known as "sphaleron" 
||1CI|| . It corresponds to a field configuration that starts at the false vacuum almost 
reaches the true vacuum 0+, and returns to 0_. It looks like a "droplet" and its explicit 
form is given by [^,11 



(t)sph{x) = 0- + 
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Sn = -cosh ^ — 
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(2.4) 
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The parameter Xc.m. reflects the translational invariance of the equation of motion, e 
gives a measure of the difference in depth between the true and the false vacua and sq is 
related to what we may call the radius of the "sphaleron" (see Figure 2). As we will see 
later, this latter parameter assumes great importance in the analysis of stability. 



It proves convenient to expand the field and and its conjugated momentum 7r{x) 
around the "sphaleron" configuration, so as to express the Hamiltonian in this particular 
basis. Collecting terms up to second order, the Hamiltonian acquires the following form 

n = E,,, + f + ^ - ^ + ^ TX-i + -r^]) (2.6) 

where {cj^} is the set of eigenvalues of the fiuctuation operator. This form of the Hamilto- 
nian clearly isolates the contributions coming from the energy of the "sphaleron" itself, the 
collection of harmonic oscilators, the translational mode and, most important, the presence 
of an inverted oscillator with frequency 1). This last contribution signals the existence of 
an unstable direction in funcional space (associated with which itself is associated with 
So 0) ready to guide the decay. Indeed, a short glance at the behavior of the energy of 
the "bubble"[| as a function of its "radius", s, shows us three possible regimes: for s < Sq, 
the "bubble" shrinks and disappears; for s > sq, the "bubble" grows without limit and, for 
s = So, we have the critical "bubble", ready to fall (see Figure 3). 

C. Time evolution and decay rate 

As we have seen in the last section, the decay of our metastable vacuum may take 
place through the nucleation, via thermal activation, of "bubbles" with radius larger than 
the critical So. In order to obtain the decay rate for this process, we must calculate the 
probability current along the unstable direction, at the saddle point. 

Assuming as our initial density matrix one that represents a set of harmonic oscillators 
of frequency Uk centered at the metastable minimum 0_, and integrating over all the degrees 
of freedom but / = — 1, one obtains, after some algebraic manipulations, the reduced initial 
density matrix (projected onto the unstable direction, Z = — 1, of functional space 

^We define here a new structure, that we call "bubble", that generalizes the "sphaleron" in the 
sense that the parameter sq is promoted to the status of a dynamical variable called s. 
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where 



pr{V-i,v'-i) = -^exp \ [a{r]\ + 7]'\) + 2-fr]_iri'^-^ 



(2.7) 



AT = JE. 



-'sph 



det{K/7r) 



1/2 



(2.8) 
(2.9) 

(2.10) 



_det{K/7T) 

and 7]i = (f)i ((/)_ - </)sph(x - xo) = //(a; - a;o) 0; ; 0o = 0, see [|]). The kernels are 
defined in the Appendix. 

The time evolution of this matrix is dictated by the reduced Hamiltonian 
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(2.11) 



through the Liouville equation. Using the most general gaussian ansatz compatible with 
unitarity for the time-evolved density matrix, one obtains for the decay rate per unit volume 
(see i) 

m 



J[saddle\ = J[i]-i]ri_i=o 



1 



fi0_i(O) A{t) ^Esph Ar(0)exp^ --02_l(o)[7^(«(O)) +7(0)] B{t) 



(2.12) 



where 



A{t) 



sinh (fit) 
"fi^ [cosh2(fit) + -^sinh2(fit)]3/2 



(2.13) 



B{t) 



1 + 



tanh^(fit) 



(2.14) 



The form of the decay rate as a function of time is given in Figure 4. To understand this 
plot, we should remember our initial hypothesis. We have defined a particular initial state 
for our system and decided to observe its "natural" evolution, i.e., an out-of-equilibrium 



evolution which does not replenish the false vacuum. Therefore, the insignificant overlap of 
the initial state with the saddle point forces r{t) to vanish at t = 0. On the other hand, 
for t > we have the decay to the true vacuum. Thus, r(t) must vanish also for t oo. 
Further discussions of these results can be found in the literature 

III. INCLUDING FERMIONS 
A. Effective action 

The motivation for the inclusion of fermions in our scheme went beyond a simple gener- 
alization of the formalism. It was based on the belief that they would introduce qualitatively 
different results for the stability of the "bubble" structures. This hope relied on the Pauli ex- 
clusion principle. As a consequence of the repulsive interaction of the fermions, we expected 
to find non-trivial metastable minima in the function E = E[s]. 

Throughout this section, we will be interested on the effects of fermions on the bosonic 
field. The natural way to proceed is to obtain, from an original fermion-boson Lagrangian, 
an effective theory for the bosons. The original Lagrangian is defined as 

C = l{d,<f>W<jy) - [V{<jy) - V{<f>,)] + ^- /i - #) (3.1) 

where /i is the bare mass of the fermions, g is the coupling constant, ipa{x) is the fermion 
field, a denotes fermion species and 02 is a local minimum of the potential 

V^(0) = y(0-0o)'(0 + 0o + y) +J0 (3.2) 

where 0o is a constant and j is an external current, responsible for the asymmetry of the 
potential even in the purely bosonic case. 

We may now integrate the generating functional 



^We may find a physical realization of this form of potential in the description of conducting 
polymers [12|. 
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_ \ fT rex, 

Z = I [Dija] Wa] exp Z / dt dx £[0, tPa, H 

Jo J-oo 



(3.3) 



over the fermions, in order to obtain an effective action for the bosons. Following the 
methods of D ashen et al. |1^, we arrive at the following expression for the effective action 

■1 



Seff[4>\ 



T roo 

dt / dx 



+ Nj2a,[<P]-J2n,a,[(P] (3.4) 



where N represent the number of fermion "species" , rij the occupation number of the discrete 
(bound states) fermionic levels and are the Floquet indices, defined such that 



i;{x,t + T) = e-^"^^(x,t) 



(3.5) 



Assuming the existence of only two symmetric bound states and the possibility of "dop- 
ing", which will characterize particular configurations defined by the occupation numbers 
n±i of the bound states, we may write the effective action as 



S, 



eff[ 



dt / dx 

"'-oo 

— dt dx 

Jo J-oo 



+ 



^(02-0^) + 6V(0-02) 



+ 



+ ^ {ai[(t>] - ai[(t>2]) - (r2+ - n_)(ai[</)] - ai[(f)2]) 



(3.6) 



i<-l 



where we have already introduced counterterms and renormalization conditions ||14||(that 
define b* through Zi(02, A, (yf, /i, j) = b*g'^). Note that n± represents the occupation number 
of particle states only, as opposed to n±i which represents particle and anti-particle states 



(n+ = n^i] = N — 



B. SPA approximation 

In order to obtain time-independent solutions, analogous to the "sphaleron", it proves 
useful to rewrite the action (|3.6| ) in terms of scattering data. Using well known results from 
inverse scattering methods |[T^-[16|, we may write the effective action as 
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s, 



effi 



[h - nlo) - b*gIo 



9 



/1--/0- 
9 



TT 



(n+ - n_)[u;o(0) - t^o(02)] 



(3.7) 



where Ui = y Jq, /i, I2 and Is have their exphcit forms shown in the Appendix. Now, 
we may extremize the action with respect to the scattering data, namely, the reflection 
coefficients r±{k) and Kq = (m|, — cu^Y^'^, 







(3.8) 



6K^ 6K^ 







(3.9) 



From ( ^.8| ) we obtain r^{k) = r_{k) = 0, ie, the potential is reflectionless. Using this, we 
may write the energy as 

2 



8 / mp 

3 It. 



N A 

TT irtp 



E = mF{-\ — \ sin + 2 ( — ) sin + ( 1 - — + ^ 1 cos^n - 



2^0 An 



In 



sin 6*0 ' 



(A^ + An) 



KmpJ \^1 + sin ^^0/ ) 
and the equation for the extrema (the gap equation) as 

TT / A \ ^ 27r / niF 
Oo + -[ tan^o ■ 



N V 9 

where we have introduced the following definitions 



sin 29q 



TT 



An 



h2)=l^ + 9'p2 = {Kl + ulYl' 



(3.10) 



(3.11) 



(3.12) 



An = n+ — n_ 



(3.13) 



ujq = mp cos ^0 



(3.14) 



A = fi-gb* -j/g 
8 



(3.15) 



Equation determines the possible values for Kq. From these, together with the 

Gel'fand-Levitan-Marchenko equation [|T^-|T6|, we may obtain the explicit form of (pbubUei^) 

K"^ 1 

4>huhhie{x) = 02 . . — T {tanh \Kq{x + So) + 5] - tanh \Kq{x - sq) + 5]} (3.16) 

gujQ smli (iiloSoj 

where 

i^o^o - tanh- (^^^) (3.17) 

tanh^^ 5°;"^"-"°|-^°"° (3.18) 
and Co is the normalization constant of the bound states. Equation (|3.16| ) may be rewritten 



as 



(pbubbUx) = 02 + 0p{tanh(^ + ^o) - tanh(^ - ^o)} (3.19) 



where 



^ = Kox + 6 (3.20) 
eo = Koso (3.21) 

(3.22) 



-"■0 



^ gujQ sinh (2i^oSo) 
Looking at ( 3.19|) we notice that, in spite of the presence of fermions, the functional 



form of the time-independent solution remained unchanged. The effect of the fermions is to 
rescale the parameters of the "bubble" . 



C. Stability of the "bubble" solution 



In order to have an asymmetrical potential with a false vacuum located at positive 02 



(see Figure 5), we must adjust the parameters such that [T| 



2g(j)e{4>e + a){(f)e - a) = - gb* - j / g = A < 



(3.23) 
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where cj) = cj) + fi/g, a = (</)q — 1/2)-^/^ and (pe is an extremum of V{(f>). 

We are now able to obtain, numerically, the solutions of the gap equation ( p.ll|) and the 



form of the energy of the "bubble" as a function of its radius. Assuming = 2, we can plot 
these results for the interesting cases An = —2, —1, (see Figures 6 and 7). 

Looking at these figures, we immediately see that our expectations are fulfilled in several 
cases. The reader will have noticed that the existence of a metastable "bubble" is associated 
with "doping": for An = —2 ("ground state") we find just the "sphaleron". 

If we remember that we have an effective theory for the bosons, it should be clear that we 
are in the same situation as at the end of section 2.B. The differences are in the rescaling of 
the parameters of 4>bubbie{x) and in the form of E = E[s]. This latter difference will certainly 
imply a difference between the lifetimes of the "sphaleron" (the "bubble" of critical radius) 
and that of the metastable "bubble" . The decay rate may be calculated in the same manner 
(see equation (|2.12| )) with E^ph replaced with Egph — Emb, Emb being the energy of the 



metastable "bubble", whenever appropriate. 

We should remark that, for temperatures which are low compared with Esph/kb, the 
metastable "bubbles" may have an appreciable lifetime. Thus, if we consider small oscilla- 
tions of these structures around the metastable minimum of E = E[s], ie, oscillations in the 
size of the "bubble" , we may be able to detect their presence through their charges 

IV. FINITE TEMPERATURE AND CHEMICAL POTENTIAL EFFECTS 

A natural way to control the "doping" mentioned in the last section consists in the 
introduction of a chemical potential, ep- This may be implemented by adding, to our 
original Lagrangian, a term like |1T7[ 



'Cch.pot. = epiploic (4.1) 



^Natural systems for this kind of experiment are doped linear polymer chains [12|. 
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The net effect of this term is to alter the form of our Floquet indices to 



(4.2) 



We may now study the case of finite temperature by doing the analytic continuation 
iT — > /3, so that 



(4.3) 



Using the results of D ashen et al. [13| for the generating functional of the fermions 



Z 



n (1+ 



^ N K/2) 

g i= — oo 



we may calculate the free energy F = —^\nZ 



F = -N{uJo - mp) - 2N [ 

Jo 



A dk 

2^ 



iu - ep - 



(3 



In 



J2 C{{ni},N) e-/5('^+i+'^-i)('^o-^F) 
.{"±i} 



and, finally, the energy due to the fermions, Ep = P ^ + F, 

Ep = ((n+) - (n_))(c.o - mp) + ^ - AN 



where we have introduced the following definitions in expression 

(a) 1'^* term: contribution from the discrete levels with 

^""'^ En.C{{n,},N)e-f^M^o-e,) 

(b) 2""^ term: contribution from the "Dirac sea" . 

(c) 3^'^ term: contribution from the continuum states with 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



A = 2/5 



du 



mp 



tan 



uo^ — nil 



(4.8) 
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where fl is the cutoff value of the energy (associated to 
the bandwidth in polymer models). 

The energy Ep may be rewritten in terms of the parameter 60 as follows 

2N 



Ep = mp { (("-+) - ("'-))(cos^o ~ 1) + 



TT 



TT 



9q cos 610+ 



where we have introduced the parameter 



AN OA 

mp dp 



(4.9) 



^ ( l^-gb* -j/9 \ ^ ^ A ^| 



(4.10) 



Ng'^ \ n + g(j)2 J Ng"^ \mp J 

which defines the form of the potential. Adding to ( [4.9| ) the bosonic contribution for the 
energy of the "bubble" we obtain, as our total energy, 



E = mp 



8 / mp 

3 It, 



V vr mpj \ TT N J 



\mp J V 1 + ^0 



(Ar+ (An)) 



4A^ (9A 



(4.11) 



^o + 7tan^o + ^ ( — j sm 2^0 = 2 



(An) 4 d d ~ 



(4.12) 



mp 8(3 J 

where we have defined (An) = (n_|_) — (n_), and the gap equation 

sin 9/9n = - 1 -^ A 

m^r d(3 

to be compared with ( p.ll| ). We are now ready to analyze the form of = E[s\ in the limits 
of high and low temperatures. 

In the high temperature limit, I3\VL — ep\ -C 1. Thus, we can approximate the values of 
A and (An) by 



Ak. (3 duo tan"^ 



u"^ — mp ^ 



p- 



(4.13) 



(An) 



in order to obtain, for the total energy. 



(3{uo - ep) 
1 - j3{ujQ - ep) 



(4.14) 
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1 — P{uJo — ep) J KmpJ \^ + sin 6*0 



'n-1 ^""-^^^ ]-AN 



- 7 sm 6'o - 



TT 

- - 6*0 ) cos + 



(4.15) 



The form of E = E[s] is shown, for various values of the chemical potential, in Figure 8. 
From this figure, we clearly see that high temperatures prevent the appearance of metastable 
"bubbles" , just yelding dissociation curves. 

For low temperatures, we have l3\mp — epl ^1. The curves for E = E[s] are shown in 
Figure 9 for the same values of e^^. In this case, temperature plays a small role, as kT is 
small compared to the other energy scales. Our results are, thus, analogous to the ones we 
encountered when we considered "doping" via a choice of occupation numbers. 

V. CONCLUSION 

The aim of this work was the study of metastable systems, with particular interest in 
the stability of "bubble" structures. The novelty, in comparison with other recent papers, 
is the presence of fermions and the exciting role they play. As we have seen, the presence 
of fermions creates a qualitatively different scenario for the evolution of the "bubbles". 
Depending almost exclusively on "doping" (relative occupation of bound fermionic states), 
which may be implemented by means of a chemical potential, we may find different regimes 
for the stability of the "bubbles". More precisely, for non-trivial "doping" we encounter 
a phenomenon of "quantum stabilization" that brings about metastable "bubbles". This 



^For ep < 0, we have the same results with anti-particles in the place of particles. 
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striking feature, that comes as a consequence of the introduction of fermions, may cause 
important changes in the physics of the models used to describe the systems mentioned in 
the Introduction. 

These results are a great stimulus for generalizing this work to higher dimensions. How- 
ever, there are some points in our scheme that require improvement: the inclusion of spinodal 
decomposition to compete with the nucleation mechanism, and a careful study of tunneling, 
as a competitor of thermal activation, together with an improvement of some approxima- 
tions, such as the use of the inverted oscilator [p|,|9|,|l8| , are some of the points to be refined 
in order to bring the formalism closer to reality. 



Acknowledgements 

E.S.F. would like to thank R.M. Cavalcanti and D. Barci for stimulating discussions. 
E.S.F. and C.A.A.C. would like to thank D. Boyanovsky for very fruitful conversations and 
collaboration. The authors acknowledge CNPq for finantial support. 



14 



f dk -ik{x-y) 


ujk cosh.{l3huJk) 


J 2n 


h sinh {I3hujk) 


f ^ p-ik(x-y) 




J 27r 


h sinh {(3hujj.) 



Appendix 

(a) Definition of tlie kernels: 
Ki{x -y) = 

K2{x -y) = 

{Ki)u, = J dxdy ^ e-'^^''-y^ Ki{k) fi{x - Xc.m.) fi'iv - Xc.m.) 

Ku' = {Ki-K2)w ; IJ'^-l 
Qi = {Ki-K2)-u; l^-l 

wliere {fi{x)} are tlie eigenfunctions of tlie fluctuation operator. 



(b) Definition of the inverse scattering integrals |14 



1 

/o ^ —V / dq 

Am J-oo 



P+(g) ^_P-{q) 



imp + q imp — q 



+ In 



/ nip — Kq 



Kmp + Kq 

^^^-\{t^j dl + 2^0+ + [(+) ^ (-)]} 



1 kdk 



2 \Jo 



m. 



\ J-oo k — 



— + 2arctan i — — 
q \ k 



[(+) 



1 r / /"oo 8 
/3 = -2 [-V y_ dq q' P^g) - ^-{K^f + [(+) ^ (-)] 

where V means principal value, = Kq and P±{q) = ln(l — |r±(g)p 
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Figure Captions: 

Figure 1: Form of the potential. 

Figure 2: The "sphaleron" solution. 

Figure 3: "Bubble" energy as a function of its radius. 

Figure 4: Decay rate as a function of time. 

Figure 5: Form of [V{(f)) - V{(f)2)]- 

Figure 6: Numerical solution of ( 3.11 ) for (a)An = 0, {h)An = —1 and (c)An = —2. 



Figure 7: Energy of the "bubble" as a function of its radius for (a) An = 0, (b) An = — 1 



and {c)An = —2. Equation (|3.17|) imposes the restriction sq > l/2mp. 



Figure 8: Energy of the "bubble" as a function of its radius in the limit of high tem- 
perature. (a) mp < ep < Q. (b) uq < ep < mp. (c) < < c^o. 

Figure 9: Energy of the "bubble" as a function of its radius in the limit of low temper- 
ature. (a) nip < ep <VL. (b) ujq < tp < mp. (c) Q < tp < ojq. 
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